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INTRODUCTION 


The  objective  of  this  work  was  to  provide  a  mathematical  model  of  the  un¬ 
winder  ribbon  capable  of  predicting  the  performance  of  actual  hardware  using 
known  parameters  as  input  data. 


This  analysis  differs  from  previous  work*  in  that  it  considers  a  moment 
equation  derived  from  the  large  deflection  theory  rather  than  the  conventional 
small  deflection  formulations.  A  geometry  change  in  the  bridge  between  the 
coiled  (unstressed)  and  the  recoiled  (stressed)  material  is  also  considered  and 
is  now  assumed  to  be  a  half -circle  rather  than  a  straight  line,* 2  As  a  result  of 
the  modifications  of  the  mathematical  model,  fall-off  (malfunction)  behavior  can 
now  be  predicted. 


UNWINDER  DEVICE 


The  basic  components  of  the  unwinder  device  are  shown  schematically  in  fig¬ 
ure  1.  The  spring  (A)  is  wrapped  around  and  fastened  at  its  inner  end  to  the 
shaft  (B).  The  outer  end  of  the  spring  is  fastened  to  the  outer  case  (C)  at 
point  D.  The  outer  case  (C)  is  fixed  to  and  rotates  with  the  projectile.  The 
axis  of  the  spring  and  the  inner  shaft  (B)  are  coincident  with  the  longitudinal 
axis  of  the  projectile.  Upon  firing,  torsional  acceleration  causes  the  spring  to 
wind  up  tightly.  After  the  torsional  acceleration  ceases,  the  centrifugal  forces 
acting  on  the  spring  will  unwind  it.  The  bending  moment  in  the  spring  bridge  and 
friction  at  the  hub  axis  inhibits  the  unwinding  process.  During  this  unwinding 
process,  the  inner  shaft  (B)  will  rotate  relative  to  the  housing.  This  motion 
can  be  used  to  close  a  switch,  rotate  a  firing  pin  in  line  with  a  detonator,  or 
initiate  some  other  arming  process. 


ANALYSIS 


Coordinates  and  Position  Vectors  (fig.  2) 


For  a  constant  angular  velocity,  w  ,  define  (xQ,  yQ,  zQ),  (x,y,z)  =  inertial 
frame  of  reference  and  body-fixed  (imbedded  in  the  hub)  coordinate  axes,  respec¬ 
tively,  so  that 


[x(t=0) ,y(t=0) ,z(t=0) ]  =  [x0,y0,z0] 


(1) 


^  "A  New  Analysis  of  the  Unwinding  Ribbon  as  a  Delayed  Arming  Device,"  W.  P. 
Dunn,  Proceedings  of  the  1982  Army  Science  Conference. 

2 

The  unwinding  ribbon  is  a  spiral-wrapped  spring  made  from  flat  metal  stock 
closely  wound.  In  the  unstressed  condition,  all  coils  of  the  spring  are 
touching. 
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Figure  1.  Delayed  arming  unwinder  device 


% 
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Figure  2.  Basic  geometry  and  physical  characterization 


» 
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The  unit  vectors  (i,j,k)  of  the  imbedded  coordinate  system  are  related  to  the 
inertial  reference  unit  vectors  (I,J,K)  by 

A  A  A 

i  =  cos(w  t-9)I  -  sin(w  t-0)J 

o  o 

A  A  A 

j  “  sin(co  t-0)I  +  cos(w  t-0)J 

o  o 

A  A 

k  -  K  (2) 

where  0  is  the  angle,  the  hub  rotates  with  respect  to  the  casing  and  is  taken  to 
be  negative;  t  is  time  and  a  is  the  angle  characterizing  the  amount  of  ribbon 
material  wrapped  onto  the  casing.  The  position  vector  to  the  point  at  which  the 
center  line  of  the  ribbon  material  just  uncoils  from  the  spool,  r,  and  the  posi¬ 
tion  vector  to  the  point  at  which  the  center  line  of  the  ribbon  material  just 
begins  to  coil  onto  the  housing,  r^  are  defined  by 

*  A  A 

r  =  -  r[cos(a-0)i  +  sin(a-0)j] 

-A  A 

r2  ”  r2 l cos (a“9 )i  +  sin(a-Q)j]  (3) 


where 


r  =  rQ-6(a-0)/2TT 

r2  =  r20~5a^2ir  (4) 

6  is  the  ribbon  thickness,  and  rQ>r20  are  initial  lengths  of  vectors  r  and 
r^,  respectively 

A 

It  =  (  V®):L  (5) 

"  dt 


For  later  use,  note  that 


di 

dt 


-(»0-o)j 


where 


d0 


The  position  vector  to  the  center  of  gravity  of  the  half -circle  bridge, 
is  defined  by 


„AB 


=  2 


2  ~  1 

+  r  +  —  k  x  -  (r  -r) 

TT  2  ,2  „ 


(6) 


4 


Substituting  equation  3  into  equation  6,  performing  the  calculation,  and  simpli¬ 
fying  produces 


rAB  =  [  T  (r2~r)cos(a-0)  -  -jj-  (r2+r)sin(a-0) 

+  [  4-  (r„-r)sin(a-0)  +—  (r„+r)cos(a-0) 

Z  Z  7T  Z 

To  obtain  a  relationship  between  a  and  0,  the  length  of  ribbon  material 
unwrapped  from  the  spool  requires  an  equal  length  wrapped  onto  the  casing,  assum¬ 
ing  the  bridge  material  dimensions  remain  constant.  Therefore, 

(a-0)  a 

/  rd(a-e)  =  /  r^da  (8) 

o  o 

Substituting  equation  4  into  equation  8,  integrating,  and  simplifying  yields 


]i 

]i 


(7) 


a 


r 

o 


2 


(9) 


Equations  of  Motion 


The  vector  quantity  ¥'  (fig.  2)  represents  the  centrifugal  force  acting  on 
the  half -circle  bridge  at~  its  center  of  mass.  Again,  assuming  the  mass  of  the 
bridge  to  be  constant, 


F 1  - 
1AB 


=  -m 


(10) 


where 


m 


AB 


-  ’5SPb5RAB 


(11) 


is  the  mass  of  the  ribbon  bridge,  b  is  the  width,  X,  is  the  density,  and  R.„  is 
the  radius  defined  by  SP 


RAB  =  l/2  <r2  +  r)  02) 

Since  there  is  only  rotational  motion,  consideration  of  the  forces  acting  on 
the  free-body  of  the  hub  and  wound  coils  (fig.  3)  yields 


* 


5 


/ 


Figure  3. 


Free-body  diagram  of  the  hub  and  uncoiled  ribbon 


6 


£  F^x  =  0  =>  F  -  N^Ccosn  -  Muslim) 

£  F  =  0  =>  F  =  N  (sinr)  +  u  cosn) 

sy  In  H 


E  -  81  ->  10  =  M  -  rF  + 


(13) 


where  5^,  5  are  reference  coordinates  parallel  to  F  and  F^,  respectively.  From 
the  first  t^o  parts  of  equation  13 


cosn  -  p  sinn 
H 


sinn  —  P  cosn 

H 


or 


Fi  "  V 

Sinn  =  (77-y)  cosn 

Yi 

2  2 

Using  sin  0  +  cos  n  =  1,  the  square  of  equation  14  produces 


(14) 


smn  = 


F  -  u  F 

1  V 


cosn  * 


F  +  V, 


y  (i  +  uh2)(f2  +  f 


T~ 2. 
1 


~H  ■  -j.  ,  Ji  l  +  nH  )(F+  F 

Substituting  the  first  part  of  equation  13  into  the  last  part  produces 


(15) 


19  =  M  -  (r  - 


u  r 
mH  H 


: — )F 

n  -n  ' 


cosn  -  y„sinn 

ci 


06) 


Substituting  equation  15  into  equation  16  then  yields 


10  -  M  +  rF  =  pHrH 


07) 


where  M  is  the  bending  moment  in  the  ribbon  (to  be  derived  later),  is  the 

normal  force  of  the  hub  on  the  hub  casing,  and  F  is  a  component  of  F'  .  The 

terra  r^  denotes  the  radius  of  the  spool  hub;  P  denotes  the  friction  coefficient 

H 

between  the  hub  and  hub  casing;  I  is  the  mass  moment  of  inertia  of  the  hub  plus 

••  Q 

wound  coils  about  the  coordinate  origin;  0  =  — ^  is  the  angular  acceleration  of 

the  hub,  and  F  is  a  scalar  component  of  which  is  perpendicular  to  F,  and 

^Ad  _L 

tangent  to  the  coils  at  A.  These  terms  (F,  F^)  are  obtained  explicitly  by  using 
the  relations 


7 


(k  x  r2//r2) 


(18) 


F  =  F' 
~AB 


F  =  F' 

1  ~AB 


l2//r2 


Return  to  equation  7  and  define  and  R2  so  that 


I'aB  =  V  +  V 

where 

R,  =  4  (r„-r)cos(ot-0)  -  —  (r o+r)sin(a-0) 

12  2  it  l 

R9  =  4-  (r„-r)sin(a-0)  +—  (r„+r)cos (a-0) 

III  Tf  2 

Taking  the  time  derivatives  of  equation  20  results  in 

R2  =  (a-0)R2  • 

R2  =  («-e)R1 

•  •  •  • 

where  r,  r~  terms  are  neglected  compared  to  a,  0  terms.  One  more  time 
of  equation  21  produces 

R1  =  -(a-0)R2  -  (a-0)2RL 
R2  =  (a-0)R1  -  (a-0)2R2 

Therefore,  using  equations  5  and  21,  equation  19  yields 


Ika  =  (“0“a>^R2i  ' 


The  time  derivative  of  equation  23  then  becomes 


r'  =  -  {[aR  +  (a)  -a)2R  ]i  +  [-aR  +  (w  -a)2R  ]j} 
~Ad  Z  O  1  I  o  z 


where,  from  equation  9, 


r  v  j.  (Se'i  i  r69'!  2 

(r20"V(ro  +  2^  -  8  hr)  . 


a  =  - 


a  =  -  9  - 


(r?n  ■  rn  50/2lT)' 


"20  o 

20  "o 


«(r,„2  -  r  2) 


0  2tr  (r20  -  rQ  -  <59/2tt) 


(19) 


(20) 


(21) 

derivative 

(22) 

(23) 

(24) 


(25) 
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Substituting  equation  24  into  equation  10,  substituting  this  result  into 
equation  18,  and  simplifying  produces 


m 


F  =  [^(r 2~r )a  "  2(r2+r ) (o^-ot) 2 ] 


F±  =  2~  [2(r2+r)a  +  Tr(r2“*r)  (mQ-a) 2] 


(26) 


Defining 


Gi  E 


G2  E 


(r20  ~  ro)(ro  +  69/2lT)  ~  8  (69/lT)2 

(r20  ‘  ro  "  69/2ir)2 

6(r202  "  ro2^ 


(27) 


2Tr(r2o  “  rQ  “  69/ 2tt) ' 


equation  25  can  be  written 


a  «  -G  0 


a  =  -  GjQ  -  G20 


(28) 


The  mass  moment  of  inertia,  I,  of  the  hub  plus  wound  ribbon  can  be  written  as 

1  “  IH  +  Y  ^  CSPb(r  (29) 


10 


where  1^  is  the  measured  mass  moment  of  inertia  of  the  hub. 


Bending  Moment  in  the  Ribbon  Bridge 


It  is  assumed  that  the  bending  moment  (produced  by  a  ribbon  element  of 
radius  r  in  the  unstressed  state,  deformed  to  a  radius  of  curvature  r2  in  the 
stressed  state)  is  equivalent  to  the  bending  moment  in  a  quarter  circle,  canti¬ 
lever  beam  going  from  a  radius  r  in  the  undeformed  state  to  a  radius  r2  in  the 
deformed  state.  This  is  represented  in  figure  4. 


Figure  4.  Large  deflection  of  quarter  circle-cantilever  beam 
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Assuming  the  beam  centerline  to  be  the  neutral  axis,  the  length  of  the  neu¬ 
tral  fiber  in  the  deformed  and  undeformed  configurations  must  be  equal.  There¬ 
fore,  from  figure  4 

irr/2  =  nr2  => 

irr 

11  "  2r2 

The  bending  moment  due  to  fiber  stress  is  defined  by 

5/2 

M  =  b  /  yady 
-5/2 

where  the  integration  is  over  the  ribbon  thickness  and  a,  the  stress  in  the 
fibers,  is  related  to  the  fiber  strains,  e,  by 

a  =  Ee  .  (32) 

Where  E  is  the  elastic  modulus  and  e,  the  strain,  is  defined  by 


(30) 


(31) 


e 


AL 

L 


(33) 


Here  L  is  the  original  fiber  length  and  AL  is  the  change  in  fiber  length  due  to 
deformation.  They  are  defined  by 


L  =  j  (r+y )  AL  =  j  (r+y)  -  n(r2+y)  (34) 

Substitute  equation  34  into  33  and  obtain  from  equation  32 

a  =  E  (l-r/r2)  (35) 

Then  substituting  equation  35  into  equation  31,  integrating,  and  simplifying 
results  in 

M  =  bEr ( 1  -£-)[rln(§±§)  -  6]  (36) 


Computational  Fora  of  the  Equations  of  Motion 


Having  the  necessary  ingredients  to  evaluate  equation  17,  the  various  terms 
will  be  redefined  for  convenient  handling  and  computational  facilitation.  Sub- 


10 


4 


stituting  equations  28  into  26  and  substituting  these  results  into  equation  17, 
solving  for  0,  and  simplifying  produces 


where 


9  '  2\  (A?  V  A7  -  4Vs  ’ 


(37) 


2  2 

A  =  a  -  rH  A 

A8  A3  ,  .  2  A5 


1  +  u 


a  a  (UHrHGl}  A 

A6  A1  ,  ,  2  A4 


H 


1  +  U 


H 


A7  "  A2  + 


2  yH  rH  G1G2  .  *2 


ry  / 

„  A,  0“  A_  =  A.  [G,  0  +  (to  +G.0)  1 

2  4  5  4L2  o  1  J 


1  + 


“H 


mAB  2  r  2 


AA  =  ^  (r?_r)  +  4  (r9+r)  ] 


(38) 


a3  =  {M - -  [ir(r2-r)G202  +  2(r2+r )  (uM^  6)2  ] } 

A2  =  {M  _  v  mABr(r2+r)(a,o'K;i®)2  "  1  mABr(r2_r)G2®2^  X 

{2I+mABr(r2-r)Gi} 

A1  =  [I  +lr(r2"r)mABGl]2 


Equation  37  is  the  defining  equation  of  motion  of  the  wound  coils  of  the  unwinder 
plus  the  hub.  The  particular  solution  is  obtained  from  the  initial  conditions 

9 ( t  =  0)  =  0  0  (t  =  0)  =  0  (39) 


Supplemental  Computations 


To  complete  the  computations,  the  values  of  rQ  and  a  Ay ,  0^^.  (the  angles 
denoting  a  completely  unwound  ribbon)  must  be  determined.  The  value  r  takes  when 
the  ribbon  is  completely  unwound  (eq  3)  is  r-^Q  and  therefore 


ri0  ro  ~  2 it  ^MAX  "  6MAX') 


(40) 


* 


11 


Also,  for  a  given  ribbon  length,  L, 

“max  “  9max 

L  =  /  rd(cx-e)  =  “  eMAX)  “  4^  (aMAX  “  0MAX) 

O 


(41) 


Solving  for  a  -  9  from  equation  40;  substituting  the  result  into  equation 
MAX  MAX 

41,  and  solving  for  rQ 


r 

o 


r  2  SL\  1/2 

^rio 


(42) 


The  maximum  wrapping  angle, 


a  ,  is  obtained  by  noting  that 


“max 

L  =  /  r^da  = 
o 


_ 6  2 

r20aMAX  “  4 it  “MAX 


Solving  equation  43  for  a  results  in 

MAX 


(43) 


MAX 


(44) 


Substituting  equations  42  and  44  into  equation  40  gives  the  radians  the  hub  ro¬ 
tates  relative  to  the  fuze  casing  as 


°MAX 


SL/ TTr 


20 


2) 


+  r 


10 


(i 


1  -  SL/  Trr 


“)} 

10 


(45) 


RESULTS  OF  COMPUTER  SIMULATION 


The  logic  and  equation  of  motion  for  the  arming  time  of  the  unwinding  ribbon 
wave  was  programmed  for  the  CDC  6600  computer  using  FORTRAN,  A  fourth  order 
Runge-Kutta  routine  was  used  to  solve  the  differential  equation  37,  The  basic 
program  and  sample  input-output  for  spring  1  is  contained  in  the  appendix. 
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A  if  jThj3 result®  of  this  analysis  compared  to  the  experimental  findines  of  T  R 

sho6  ShOWn  ln  f,igures  5  through  11.  The  analysis  illustrated2 in  figures 
and  11  shows  more  resistance  to  motion  than  indicated  by  experiment.  ^par¬ 
ticular,  where  none  of  the  ribbons  completely  unwind  (fig.  5)  the  analyst  is 

X7oatlVIM  th8t  ^  “I1S  are  PredlCted  to  the  con  then act" 

ally  do.  This  error  appears  to  be  less  as  the  angular  speed  increases-  therP- 

^  °VerWhelmS  b°th  the  bendi^  movement  in  the 


A  parametric  study  was  made  to  determine  the  pffppf  . 

quantities  on  the  results  of  the  analysis.  The  effects  of  varying  u  r°U/r  lnpUt 

o’  10  20’ 

5,  and  yH  are  shown  in  figures  12  through  15.  It  appears  that  the  arming  time 

the  raSoenr  T/ K1S  “°8t,  sfns1itive  to  the  ratio  r10/r20  and  therefore 
ribbon  r10/r20  could  be  *  critical  factor  in  the  design  20°f  the  unwinder 


CONCLUSIONS  AND  RECOMMENDATIONS 


ful  design  CooTYorThe"  *  P°tentiill>'  po”“- 


As  adjustment  of  the  parameters  a)  r  /r  &  ,,  n n  , 

conducted  ana,lytic  and  experimental  results0,  is  ’recommended  that^new  tests *11 
fric^on  coefficients  °f  the  actual  dimensions,  angular  speeds,  and 

T^^fThetite^T.l^  ^  th* 


~AB  dt  ^mAB  ^AB  ^ 


'  “AB  ^AB  mAB  ^AB 


nnt-  K1118’  h°^evet»  would  involve  considerable  mathematical  complexity  and  should 
not  be  considered  unless  the  comparison  of  theoretical  and  new  experimental 

coned VETh^  “f8SMr  ^  “»»<*  of  the  assumption  the t The  nbW 
coiled  on  the  hub  is  in  a  stress-free  state  should  be  verified. 


Study  of  Wrapped  Springs  for  Application  to  a  Delayed  Arming  Device  "  T.  B 
lfnend,  Summary  Report  ER-1404,  Aircraft  Armaments,  Inc.,  1958. 
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Figure  5. 


Spring  1,  relative  turns  versus  time 
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1111m 


Figure  6. 


Spring  2,  relative  turns  versus  time 
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Figure  7.  Spring  3,  relative  turns  versus  time 
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Figure  8.  Spring  4,  relative  turns  versus  time 
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Figure  9.  Spring  5,  relative  turns  versus  time 
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Figure  10.  Spring  6,  relative  turns  versus  time 
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Figure  11 


Spring  7 


relative  turns  versus  time 


* 


» 


iboo  r/s. 


^c.% 


S  --C.60H 


^.=ew  v/k 
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Figure  12.  Spring  1,  effect  of  spin-rate  (co  )  on  arming 
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Figure  13.  Spring  1,  effect  of  r  0/r2Q  on  arminS  and  armin8  time 
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Figure  14.  Spring  1,  effect  of  ribbon  thickness  on  arming 
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Figure  15.  Spring  1,  effect  of  hub  friction  (y  )  on  arming 

H. 
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APPENDIX 

COMPUTER  PROGRAM  AND  SAMPLE  INPUT-OUTPUT  DATA 
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SUBROUTINE  OUTP( T I  ME , THET A , DTHET A , IHLF , NDIM , PRMT  ) 

DIMENSION  THET A( 2 ) , DTHETA ( 2 ) ,PRMT(5) 

COMMON  DELTA, B, RIO ,R20,R0, E , RHOSP , RHOSH , M , L.OMO , RC, ALPHA, PI , I 
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TIME*  .850  TURNS*  22.59  THE 


TIME*  .900  TURNS*  22.65  THETA( 1 ) =-1 42. 31  M*  .9077E-01  THDOT*  .466774E+02  THETA( 2) *-.581 903E+01  ALPHA*  68.25 
TIME*  .950  TURNS*  22.69  THETA( 1 ) =-1 42. 55  M*  .9084E-01  THDOT*  .314048E+02  THETA(2)*-.389t40E+01  ALPHA*  68.35 
TIME*1 .000  TURNS*  22.71  THETA( 1 ) =-142. 71  M*  .9089E-01  THDOT*  .21 0457E+02  THETA(2)«-.259703E+01  ALPHA*  68.41 
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